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Abstract — Under  the  National  Hurricane  Program,  MIT  Lin¬ 
coln  Laboratory  is  developing  HURREVAC  extended  (HVX), 
which  is  replacing  the  legacy  evacuation  decision-support  tool 
HURREVAC.  The  platform  is  designed  to  enable  emergency 
managers  (EMs)  make  accurate  and  timely  evacuation  decisions 
for  their  community  by  providing  planning,  situation  awareness, 
data  analytics,  and  training  tools  in  a  Web  browser.  HVX 
provides  the  novel  capability  to  simulate  storm  tracks  in  a  Web 
application  with  minimal  user  input  or  meteorological  expertise. 
This  capability  will  allow  EMs  to  train  more  frequently  and 
inexpensively,  simulate  storms  that  EMs  have  not  already  seen, 
and  enable  serious  games  and  training  analytics  to  improve  EMs 
preparedness.  The  simulator  leverages  the  HURDAT2  database, 
which  contains  information  on  more  than  2,800  tropical  storms, 
to  generate  plausible  storms.  The  intended  purpose  and  the  requi¬ 
red  storm  characteristics  drove  simulator  design  and  technique 
selection.  Storm  movement  between  waypoints  is  simulated  by 
a  Markov  model  and  kernel  density  estimation  (KDE)  with  a 
gamma-von  Mises  product  kernel.  KDE  with  a  beta  kernel  gene¬ 
rates  maximum  sustained  winds,  and  linear  regression  simulates 
minimum  central  pressure.  Maximum  significant  wind  extents 
are  simulated  by  Poisson  regression  and  temporal  filtering.  The 
un-optimized  MATLAB  code  runs  in  less  than  a  minute  and  is 
integrated  into  HVX  as  a  Java  package. 

I.  Introduction 

Hurricanes  and  typhoons —  tropical  cyclones  with  maxi¬ 
mum  sustained  winds  of  64  knots  (74  mph)  or  greater  [1] — 
can  cause  many  injuries  and  fatalities  as  well  as  millions  or  bil¬ 
lions  of  dollars  in  damages  [2].  Tropical  cyclones  are  usually 
categorized  by  their  maximum  sustained  winds  according  to 
Table  I,  with  Category  3  and  higher  storms  considered  “major.” 
Storm  surge,  inland  flooding,  and  extreme  winds  cause  the  vast 
majority  of  fatalities  and  destruction  [2]. 

Emergency  managers  (EMs)  are  responsible  for  making 
hurricane-related  evacuation  recommendations  to  higher  of¬ 
ficials.  If  ordered,  an  evacuation  must  be  completed  before 
the  onset  of  tropical  storm-force  (34-knot)  winds.  An  EM 
must  coalesce  information  about  the  community’s  hurricane 
evacuation  plans,  National  Hurricane  Center  (NHC)  weather 
forecasts,  and  information  from  local  National  Weather  Ser¬ 
vice  (NWS)  offices;  coordinate  with  emergency  services,  neig- 
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TABLE  I 

Tropical  Cyclone  Strengths 


Abbreviation 

Strength 

Maximum  Sustained 

Winds  [knots] 

TD 

Tropical  Depression 

[0,  33] 

TS 

Tropical  Storm 

[34,  49] 

STS 

Strong  Tropical  Storm 

[50,  63] 

HU1 

Category  1  Hurricane 

[64,  82] 

HU2 

Category  2  Hurricane 

[83,  95] 

HU3 

Category  3  Hurricane 

[96,  112] 

HU4 

Category  4  Hurricane 

[113,  136] 

HU5 

Category  5  Hurricane 

[137,  oo*) 

*The  Saffir-Simpson  Hurricane  Wind  Scale  (SSHWS)  [3]  does 
not  specify  an  upper  limit  on  the  maximum  sustained  winds  of 
a  Category  5  hurricane,  but  the  simulator  in  this  paper  imposes  an 
upper  limit  of  180  knots  (207  mph). 

hboring  governments,  and  relief  agencies;  and  communicate 
with  the  public. 

A.  National  Hurricane  Program 

Following  Hurricane  Katrina  in  2006,  the  U.S.  Government 
mandated  the  DHS  S&T/FEMA  National  Hurricane  Program 
(NHP).  In  2014,  NHP  launched  a  Technology  Modernization 
(TM)  effort,  which  identified  technology  gaps,  explored  new 
solutions,  established  a  Working  Group  (WG),  and  started  new 
technology  development  [4]. 

For  many  parts  of  the  continental  U.S.,  the  hurricane  return 
period  exceeds  eight  years,  and  the  major-hurricane  return 
period  exceeds  three  decades  [2].  The  last  storm  to  make 
landfall  as  a  major  hurricane  in  the  continental  U.S.  was  Wilma 
in  2005.  Consequently,  an  EM  might  be  familiar  with  the  few 
historical  storms  that  affected  his  or  her  community  but  have 
limited  real-life  evacuation  experience. 

Limited  travel  budgets,  training  sessions,  and  availability  of 
NHC  experts  mean  that  less  than  10%  of  EMs  receive  training 
each  year.  Sessions  last  1-5  days,  focus  on  hurricane  forecasts 
and  other  NHC  products,  and  cover  only  one  or  two  scenarios. 
As  a  result,  an  EM’s  training  can  become  stale,  and  there  is 
little  data  available  to  assess  the  efficacy  of  training  or  EM 
preparedness. 

The  EMs’  legacy  decision- support  tool,  HURREVAC,  is  a 
Windows  application  with  limited  decision-support  analytics 
and  training  capabilities,  and  it  is  not  Web-accessible.  Its 


simulator  has  a  tedious,  manual  entry  process  that  requires 
technical  knowledge  of  tropical  cyclones,  so  it  is  rarely  used. 

The  NHP  TM  effort  recommended  replacing  HURREVAC 
with  an  integrated,  Web-enabled  decision- support  platform  that 
provides  planning,  situational  awareness,  data  analytics,  and 
training  tools  to  enable  EMs  to  make  accurate,  timely  evacu¬ 
ation  decisions.  MIT  Lincoln  Laboratory  began  development 
of  the  new  platform,  HURREVAC  extended  (HVX),  in  2015. 

B.  Simulator  Motivation,  Requirements,  and  Usage 

One  of  HVX’s  goals  is  to  improve  EM  preparedness  by 
enabling  frequent,  inexpensive  training  on  the  EM’s  own  com¬ 
puter,  offering  rapid-play  serious  games  with  never-before- 
seen  storms,  and  providing  training  analytics  for  preparedness 
exercises.  These  desired  properties  require  a  storm  track  si¬ 
mulator,  which  plays  a  small  but  important  role  in  HVX  and 
is  the  subject  of  this  paper. 

The  simulator  does  not  need  to  provide  the  physical  accu¬ 
racy  or  fine  detail  of  NHC  supercomputer  simulations  and 
forecasting  models.  Rather,  it  needs  to  produce  plausible  storm 
behavior;  create  tracks  that  EMs  have  not  seen  before;  accept 
a  few  simple  user  inputs  without  requiring  technical  tropical 
cyclone  expertise;  and  run  in  the  HVX  Web  application  in 
under  a  minute. 

To  use  the  simulator,  a  user  clicks  on  a  regional  map  to  enter 
a  few  latitude-longitude  waypoints  that  outline  the  desired 
storm  track.  For  each  waypoint,  the  user  enters  the  desired 
strength  from  Table  I  and  the  number  of  time  steps  until  the 
next  waypoint.  Figure  la  shows  an  example  of  the  user  inputs. 
Figure  lb  shows  the  output  of  the  simulator,  which  is  covered 
in  Section  III. 

II.  Data,  Related  Work,  and  Approach 

This  section  describes  the  primary  data  source,  some  related 
work,  and  the  simulation  approach. 

A.  HURDAT2  Storm  Track  Data 

The  simulator  uses  the  HURDAT2  database  [5],  which  con¬ 
tains  the  best  track  data  from  the  NHC’s  post- storm  analysis 
of  actual  tropical  cyclones.  Tables  are  provided  for  the  North 
Atlantic  Basin  and  the  Northeast  and  North  Central  Pacific 
Basins.  The  former  (5.9  MB)  contains  data  on  about  1,800 
storms  from  1851-2015,  and  the  latter  (3.2  MB)  contains  data 
on  about  1,000  storms  from  1949-2015. 

For  each  storm,  the  table  provides  observations  at  synoptic 
time  steps  (0000,  0600,  1200,  and  1800  UTC)  and  may  include 
times  of  events  like  landfall.  Each  observation  provides  the 
latitude-longitude  location  of  the  storm  center,  the  strength , 
and  the  maximum  sustained  winds  in  5 -knot  increments.  Since 
the  1980s,  data  for  minimum  central  pressure  in  integer 
millibars  became  available,  and  since  2004,  twelve  wind  radii 
describing  the  extent  of  hurricane-force  (64+  knot),  strong 
tropical  storm-force  (50-64  knot),  and  tropical  storm-force 
(34-50  knot)  winds  at  5-nautical  mile  increments  in  four 
quadrants  became  available. 


B.  Related  Work 

Some  aspects  of  the  simulator  are  similar  to  work  by  Rumpf 
et  al.  [6]  [7]  and  Nakamura  et  al.  [8].  Rumpf  et  al.  used  a 
first-order  Markov  model  and  nonlinear  regression  and  were 
concerned  with  estimating  risk  and  losses  for  the  insurance  in¬ 
dustry.  Nakamura  et  al.  reported  longer-term  trends  in  tropical 
cyclone  behavior.  To  capture  such  trends,  they  concatenated 
randomly  sampled  track  segments  from  HURDAT2,  with  the 
sampling  probability  based  on  the  similarity  between  historical 
storms  and  the  partially-formed  simulated  storm. 

The  present  simulator  need  only  generate  plausible  storms 
but  must  have  a  modest  computational  footprint  and  execute 
quickly.  The  presence  of  longer-term  trends  led  us  to  choose 
a  hierarchical  approach,  but  we  also  wanted  the  ability  to 
generate  never-before-seen  behavior. 

C.  Notation  and  Approach 

The  simulator  generates  results  similar  to  HURDAT2.  Let 
L  be  the  total  storm  duration.  For  each  synoptic  time  step 
t  =  1,  2,  . . . ,  L,  the  simulator  produces  a  latitude-longitude 
location  £(t)\  forward  movement  u(t)  =  (p(£),  #(£)),  where 
p{t)  is  forward  speed  in  knots,  and  6{t)  is  forward  heading 
in  degrees;  atmospherics  (strength  s(t ),  winds  w(t),  pressure 
6(f));  and  a  12-element  wind  radii  vector  r{t). 

Let  M  denote  the  number  of  waypoints;  wm  the  rath 
waypoint  for  ra  =  1,  . . . ,  M;  sm  and  tm  the  strength 
and  time  step  at  wm,  respectively;  and  Lm  the  duration  or 
number  of  time  steps  until  the  next  waypoint.  Lm  does  not 
include  the  current  waypoint  wm,  and  the  terminal  waypoint 
wm  has  Lm  =  0.  Hence,  L  =  J2m= i  Lm  +  1-  A  stage 
Sm  =  {wm,  rnm+i,  sm,  sm+i,  Lm}  consists  of  a  pair  of 
successive  waypoints,  their  strengths,  and  the  stage’s  duration. 
Hence,  a  storm  has  M  waypoints  and  M  —  1  stages. 

The  indicator  /birth  (f)  equals  one  if  t  =  1  and  zero  other¬ 
wise.  For  Wm  or  i(t ),  the  indicator  or  /iand(/(f)) 

equals  unity  if  the  location  lies  over  land  and  zero  otherwise. 

We  adopted  a  hierarchical,  data-driven  approach.  From  top 
to  bottom,  the  levels  are  Season,  Storm,  Stage,  and  Time- 
step.  The  higher  levels  are  mainly  conceptual,  as  development 
focused  on  the  Time-step  level.  The  Season  level  models 
broad,  seasonal  characteristics  such  as  a  normal,  La  Nina,  or  El 
Nino  weather  pattern,  and  the  number  of  tropical  cyclones  of 
each  strength.  For  an  individual  storm,  the  Storm  level  models 
the  number  of  waypoints  M,  birth  and  death  locations  W\  and 
wmi  and  total  number  of  steps  L.  The  Stage  level  models 
the  remaining  waypoints  w 2,  ...,  wm-u  and  the  strengths 
Sm  and  durations  Lm,  m  =  1,  ...,  M  —  1.  For  each  stage 
Sm ,  the  Time- step  level  simulates  the  forward  movement, 
atmospherics,  and  wind  radii  to  produce  £(t),  s(t ),  w(t),  b(t ), 
and  r(t)  for  each  time  step  in  the  stage. 

The  simulator  was  developed  in  MATLAB  and  leverages 
HURDAT2  to  construct  data-driven  statistical  models  that  can 
generate  plausible  yet  never-before-seen  storm  behaviors.  For 
a  particular  storm  component,  it  selects  a  set  X  =  {xi,  ..., 
X]\f}  of  representative  HURDAT2  data,  assumes  the  samples 


(a)  Example  of  three  user-entered  waypoints  with  desi-  (b)  Example  of  simulated  storm  track 

red  strengths  and  durations 

Fig.  1.  Example  of  simulator  input  and  output. 


Xi  are  distributed  f{x\X)9  builds  a  fitted  model  f(x;X)9  and 
then  draws  a  simulated  realization  x  from  f(x;X). 

The  approach  can  model  dependencies  that  extend  beyond 
a  single  time  step,  and  it  can  exploit  non-causal  relations¬ 
hips.  The  models  can  model  relationships  explicitly,  such  as 
using  a  conditional  distribution  like  f(b(t)\w(t))9  or  implicitly 
through  the  data- selection  criteria  used  to  form  A,  indicated 
by  X(w(t)).  For  example,  a  Stage  can  be  specified  to  start  as 
a  TS,  end  as  a  HU2,  and  endure  for  10  steps  (60  hours),  and 
this  information  is  exploited  by  selecting  data  on  storms  that 
strengthened  in  a  similar  fashion  over  a  similar  time  period. 

1 )  Kernel  Density  Estimation:  The  simulator  uses  both  ker¬ 
nel  density  estimation  (KDE),  a  non-parametric  method,  and 
parametric  regression  methods.  Given  X,  KDE  approximates 
/( x;X)  as  f(x;X)  =  {l/N)YJiK{x,xi)  where  K(x,X{ ) 
is  a  kernel  function ,  typically  a  unimodal  probability  density 
function  (PDF)  with  its  mode  at  X{.  As  a  result,  f(x;X)  will 
be  a  valid  PDF,  too.  A  simulated  realization  x  is  generated  by 
drawing  a  sample  uniformly  from  X ,  and  then  drawing  x 
from  a  random  number  generator  distributed  K(x,xf). 

2)  Regression:  Regression  methods  assume  that  f{x\X) 
has  a  parametric  form  f{x\(x,X)9  such  as  a  normal  distribu¬ 
tion  with  a  denoting  the  mean  and  variance.  They  typically 
compute  the  maximum  likelihood  estimate  a  to  produce  the 
fitted  model  f(x;X)  =  f(x\a.9X).  Simulation  is  accomplis¬ 
hed  by  drawing  x  according  to  f(x\OL,X). 

III.  Simulator  Component  Details 

Figure  2  shows  a  block  diagram  of  the  simulator  com¬ 
ponents,  described  in  this  section.  For  each  stage,  storm 
movement  and  atmospherics  are  simulated  separately.  Next, 
the  stage  results  are  connected  into  a  single  track.  Finally,  the 
wind  radii  are  simulated  to  yield  the  simulated  storm  track. 


waypoints 


Fig.  2.  Block  diagram  of  simulator. 


A.  Movement 

The  movement  simulator  combines  KDE  with  a  first- 
order  Markov  model,  so  f(y{t)\£(t),v(t  —  1),...  ,;X)  = 
f(is(t)l£(t),is(t  —  1);X),  which  is  modeled  implicitly  as 
f(v(t);X(£(t),Ihmd(t(t)),v(t  -  1), /birth (i)))  with  KDE.  X 
is  formed  from  data  for  locations  similar  to  £{t)  including 
over-land  or  -water  status,  forward  movement  similar  to 
v(t  —  1),  and  similar  storm  birth  or  continuation  states.  The 
set  A  =  {v\9  ...,  izjv},  where  each  Vi  =  (pi,6i)  is  a 
historical  forward-movement  pair.  Then  X(£(t),  v{t  — 

1),  «/(*)),  /birth m  =  (1  /N)  J2i  K(u(t),  Vi). 


Fig.  3.  Contour  plot  of  example  movement  distribution  with  gamma- von 
Mises  KDE  from  54  historical  data  samples.  Forward  speed  corresponds  to 
radius,  and  forward  heading  corresponds  to  angle. 


Fig.  4.  Example  of  simulated  movement  for  first  stage. 


Since  forward  speed  is  nonnegative  and  forward  heading 
is  directional,  we  adopted  a  product  kernel  = 

Kg(p(t),  Pi) Kvm(6(t),  Qi ),  where  Kg  is  a  gamma  kernel  with 
2-knot  bandwidth  [9],  and  Kvm  is  a  von  Mises  kernel  with 
15 -degree  bandwidth.  Kg  is  based  on  the  gamma  distribution, 
whose  PDF  has  positive  support  and  is  unimodal.  Kvm  is  de¬ 
rived  from  the  von  Mises  distribution,  whose  PDF  is  circularly 
periodic,  unimodal  over  a  period,  and  described  by  its  mode 
(i  and  concentration  k  >  0.  As  k  0,  the  PDF  approaches 
a  uniform  distribution  over  all  directions,  and  for  n  >  2,  the 
PDF  approximates  a  Normal(/i,  ft-1/2)  wrapped  around  the 
unit  circle  [10]  [11]  [12].  A  contour  plot  of  an  example  fitted 
distribution  appears  in  Figure  3. 

This  formulation  is  natural,  intuitive,  and  consistent  with  the 
idea  of  forward  movement,  and  it  avoids  difficulties  such  as  the 
nonlinear  mapping  between  Cartesian  and  polar  coordinates 
and  the  possibility  of  “negative”  speed  that  could  arise  with  a 
2-D  normal  kernel  often  used  for  multi-dimensional  KDE. 

By  itself,  this  model  induces  a  random  walk  beginning  at 
t(tm)  =  so  the  storm  is  very  unlikely  to  be  located 

at  wm+ 1  when  t  =  £m+i.  To  fix  this  problem,  we  con- 
vexly  blend  the  KDE  movement  with  rhumb-line  movement 
straight  to  wm+i.  As  t  increments  from  tm  to  tm+ %  —  1,  we 
first  compute  the  forward  movement  z/(f)  needed  to  move 
straight  from  i(t)  to  wm+i  over  the  stage’s  remaining  time 
steps.  We  then  compute  the  blended  location  £{t  +  1)  = 
5-1  (v(t)  +  v{t))  +  (1  -  v(t))  g(£(t)  + 1/ (£))),  where  g 

and  g~x  denote  conversion  to  and  from  Earth-centered,  Earth- 
fixed  coordinates,  respectively.  The  weight  decreases  linearly 
from  rj(tm)  =  1  to  rj(tm+i  —  1)  =  0,  so  the  straight  path  gains 
influence  as  t  approaches  £m+i,  and  the  track  is  guaranteed  to 
satisfy  £(tm+ 1)  =  Figure  4  shows  an  example  of  the 

simulated  movement  for  a  single  stage. 

B.  Atmospherics 

1)  Strength :  To  simulate  the  strength  s(t)  at  each  time  step 
in  stage  <Sm,  the  simulator  maps  the  strengths  in  Table  I  to  the 


Max.  Sustained  Winds  [kt] 

Fig.  5.  Example  of  KDE  for  maximum  sustained  winds  given  a  strength  of 
HU1  and  non-birth  time  step. 

ordered  set  {1,  2,  ...,  8}  and  simply  uses  a  noisy  logistic 
function  from  to  sm+i  over  duration  Lm.  This  is  the  only 
simulator  component  that  is  not  data-driven. 

2)  Winds:  For  a  given  strength  s(t),  the  winds  must 
lie  in  a  particular  interval  (see  Table  I).  The  simula¬ 
tor  uses  KDE  and  models  f(w(t)\s(t);X)  implicitly  as 

X(s(t),  /birth (t))  =  X  = 

{wi, . . . ,  wn}  is  formed  from  wind  data  matching  the  strength 
s(t)  and  birth  or  continuation  status.  is  a  beta  kernel  with 
5-knot  bandwidth  [13],  which  remaps  the  beta  PDF  from  [0,1] 
to  an  arbitrary  closed  interval.  Storm  birth  must  occur  as  a  TD, 
TS,  or  STS,  while  continuation  has  eight  possible  strengths, 
yielding  11  models.  An  example  is  shown  in  Figure  5. 

3)  Pressure:  To  simulate  pressure  b(t )  for  a  stage 
Sm ,  we  assume  a  simple  linear  relation  between  winds 
and  pressure  [14].  Standard  linear  regression  models 
f(b(t)\w(t);  X(Sm,  4mdOm),  /iandOm+i)))  as  a  normal  dis¬ 
tribution  whose  mean  is  a  linear  function  of  w(t).  The  set 
A  =  {(wi,bi),  ...,  (rc/v,fr/v)}  is  formed  from  HURDAT2 
storm  segments  that  are  similar  to  Sm  in  terms  of  their  starting 
and  ending  strengths  and  sm+i,  duration  Lm,  and  whether 
the  waypoints  are  over  land  or  water.  In  this  way,  the  simulator 
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Fig.  6.  Example  of  linear  regression  for  minimum  central  pressure  given 
maximum  sustained  winds. 
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Fig.  8.  Example  of  connected  simulated  stages.  Wind  radii  have  not  been 
simulated  yet. 


Fig.  7.  Example  of  simulated  atmospherics  for  first  stage. 

exploits  longer-term,  non-causal  dependencies.  Data  selection 
and  linear  regression  are  conducted  online.  Figure  6  displays 
an  example  fit,  and  Figure  7  shows  the  simulated  strength, 
winds,  and  pressure  for  a  single  stage. 

C.  Stage  Connection 

The  preceding  steps  simulate  location  &{t)  and  atmospherics 
s(t),  w(t ),  b(t)  for  all  stages  and  time  steps.  The  stages 
are  concatenated  to  form  a  track  lacking  only  wind  radii;  an 
example  appears  in  Figure  8. 

D.  Wind  Radii 

For  each  wind  force  cc  in  the  ordered  set  {HU,  STS,TS} 
and  each  quadrant  q  G  {NE,  SE,  SW,  NW}  (see  Figure  9),  the 
wind  radius  is  positive  and  can  depend  upon  multiple 

factors,  such  as  w(t),  b(t ),  and  any  inner  radii  ru^q{t), 
uj'  <  uj.  The  simulator  uses  Poisson  regression  [15]  to  model 
f(ru,q(t)\w(t),b(t),ru,,q(t),  /land (£(t));X(u,q,s(t)).  Hence, 
ruj,q{t)  is  assumed  to  be  distributed  Poisson  with  parameter 
A  =  otTx ,  where  a  is  a  weight  vector,  and  x  is  a  vector  of 
w,  b ,  any  inner  radii  v)?,  and  over-land  or  -water  status.  X 
was  selected  from  data  matching  uj,  q ,  and  strength  s(t),  and 


Fig.  9.  Wind  radii  in  four  quadrants. 

a  total  of  72  different  models  were  fitted  offline.  Figure  10 
shows  example  regression  models  for  all  four  quadrants;  the 
models  have  the  circular  asymmetry  typical  of  real  hurricanes. 

After  wind  radii  have  been  simulated  for  the  entire  track, 
median  filtering,  attenuation  over  land,  and  smoothing  clean 
up  the  wind  radii  and  complete  the  simulation.  Figure  11 
displays  an  example,  and  the  final  simulated  storm  is  shown 
in  Figure  lb. 

IV.  Summary 

The  HVX  storm  track  simulator  produces  plausible,  data- 
driven  storm  tracks  to  support  EM  training,  serious  games, 
and  training  analytics.  It  accepts  a  few  simple  inputs  and 
does  not  require  detailed  knowledge  of  tropical  cyclones.  The 
MATLAB  simulator  has  been  integrated  into  HVX  as  a  Java 
package,  and  it  runs  in  under  a  minute. 

The  mild  goal  of  plausibility  requires  less  stringent  vali¬ 
dation  than  estimating  hurricane  insurance  risk  or  planning 
levee  sites.  Intermediate  versions  of  the  simulator  have  been 
positively  received  by  the  community  at  the  2015  Interagency 
Coordinating  Committee  on  Hurricanes  meeting,  the  2016  Na¬ 
tional  Hurricane  Conference,  and  recent  NHP  WG  meetings. 


SW  Quadrant  SE  Quadrant 


Min  Central  Pressure  [mbar]  Min  Central  Pressure  [mbar] 


Fig.  10.  Examples  of  Poisson  regression  for  hurricane-force  wind  radii  for 
all  four  quadrants  given  105-knot  winds,  over-water  status,  and  Category  3 
hurricane  strength.  Hollow  circles  indicate  sample  data  for  other  wind  speeds 
(100-  and  110-knot)  or  over-land  status.  Solid  circles  indicate  data  matching 
105 -knot  winds  and  over- water  status. 


change  in  forward  heading  or  rapid  strength  could  be  injected 
into  an  otherwise  ordinary  simulated  storm. 
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Fig.  11.  Example  of  simulated  wind  radii.  Other  information  (strength,  winds, 
pressure)  has  been  suppressed. 


Current  work  is  using  simulated  tracks  to  create  forecasts, 
text  advisories,  and  discussions,  and  to  reference  products 
such  as  precomputed  NHC  storm  surge  tables.  Possible  further 
technical  work  includes  statistical  validation,  selection  of  KDE 
bandwidths  by  cross-validation,  and  implementation  of  other 
hierarchy  levels  so  the  user  can  just  choose  the  landfall  loca¬ 
tion.  Finally,  the  modeling  methods  can  provide  the  probability 
of  different  storm  behaviors,  so  surprising  events  like  a  sudden 


